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Abstract 



We calculate the combination 2aQ°' — SaJ,^' (the Olsson sum rule) and 
the scattering lengths and effective ranges ai, ttj^' and &i, b'^P dispersively (with 
the Proissart-Gribov representation) using, at low energy, the phase shifts for 
TTTT scattering obtained by Colangelo, Gasser and Leutwyler (CGL) from the Roy 
equations and chiral perturbation theory, plus experiment and Regge behaviour 
at high energy, or directly, using the GGL parameters for as and bs. We find 
mismatch, both among the GGL phases themselves and with the results obtained 
from the pion form factor. This reaches the level of several (2 to 5) standard 
deviations, and is essentially independent of the details of the intermediate en- 
ergy region (0.82 < E < 1.42 GeV) and, in some cases, of the high energy 
behaviour assumed. We discuss possible reasons for this mismatch, in particular 
in connection with an alternate set of phase shifts. 
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1. Introduction 



In two remarkable recent papers, Ananthanarayaii, Colaiigclo, Gasser and Leutwylert^I and Colangelo, 
Gasser and Leutwyler^l (to be referred to as, respectively, ACGL and CGL) have used experimental 
information, analyticity and unitarity (in the form of the Roy equations'^! ) and, in CGL, chiral calculations 
to two loops, to construct what is presented as a very precise tttt scattering amplitude at low energy, 
E = §1/2 < 0.8 GeV. 

There is little doubt that the small errors claimed by CGL, at the level of very few percent, follow 
from the Roy-chiral analysis, plus chiral perturbation theory(with the assiimption of negligible higher 
order corrections), given the input scattering amplitude at high energy, say, for E > 1.42 GeV. What is 
however not so clear is that the input selected by ACGL is unique, not even that it is the more physically 
acceptable one. The question then remains, what is the effect of changing this high energy input in the 
low energy tttt ainplittide. 

In the present paper we address ourselves to the matter of the consistency of the CGL S matrix. To 
be precise, we evaluate the following quantities: the combination of SO, S2 scattering lengths 2a\^^ — ba'^^P 
{Olsson sum rule): the scattering length ai and effective range^ 6i in the P wave; and the scattering 
lengths and effective ranges Oj^^ 62^^ J = 0, 2 for the DO, D2 waves. For the ai, 61, 03^^ 63^^ we use 
the Proissart-Gribov representation.^ This presents two advantages. First of all, it was not verified in 
ACGL or CGL; therefore, it provides a novel test of the CGL phase shifts. Secondly, for ai, bi and, to 
a lesser extent, for the 03^^ the Froissart-Gribov representation is sensitive to the high energy scattering 
amplitude, precisely one of the features we want to probe. We then compare what we find with the values 
for as and 6s given by CGL themselves. For 0,1 . 61 we also compare the CGL evaluations with the results 
of a direct fit of the P wave to the pion form factor, this last a fully independent test. 

The result of our calculations (Sect. 4 here) is that the solution of CGL is not consistent with 
the results from the fit to the pion form factor or with itself (if assuming a reasonable high energy Regge 
behaviour) and the mismatch occurs essentially independently of the details of the intermediate energy 
(0.82 < < 1.42 GeV) phase shifts we use, provided they fit experiment (Subsect. 4.4), and, in some 
cases, also of assumptions on the high energy (s^/^ > 1.42 GeV) behaviour. For some of the quantities 
discussed above the disagreement reaches several (up to 5) standard deviations. For the a^^ the more 
striking discrepancy occurs for the combination ao+ — — 02^''] • This is because it corresponds to a 

combination of only isospin 1, 2 in the s, u channels so the Froissat-Gribov integral is very accurate since te 
SO wave, large and the one less well known, does not contribute. The chiral perturbation theory calculation 
for this quantity has also small errors since (to one loop) it only depends on one chiral lagrangian constant, 
I2, see below. 

The mismatch is much less severe (below the 2 a level) for the ACGL results, the main reason 
being that their errors are at least three times as large as the CGL ones. We discuss in Sect. 5 the reasons 
for the CGL mismatch, which may be due to the use by CGL of an irrealistic high energy part of the 
scattering amplitude, which distorts their low energy (a-^/-^ < 0.82 GeV) phase shifts beyond the very small 
errors implied by their assumption of negligible higher chiral perturbative corrections. 

Apart from these two sections, we present in Sect. 2 the Roy equations, in Sect. 3 the scattering 
amplitude we will use (including in particular a detailed discussion of the high energy pieces) and finish 
the paper with summary and conclusions in Sect. 6. 



^ Actually, 61 is not the effective range, though it is related to it. We use the definitions of ACGL and CGL 
for the a and 6, except that we take the dimensions of o; to be M^^'^^'^^K Here is the charged pion mass, 
M„ ~ 139.57 MeV. 

^ The method of the Proissart-Gribov representation to calculate scattering lengths and effective rages was intro- 
duced in refs. 4,5. It is also discussed in some detail in ref. 6. 
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2. The Roy equations 
2.1. Dispersion relations 

The analyticity properties of the nn scattering amphtude, F{s,t), imply that we can write a Cauchy 
representation for it, fixing t and allowing s to be complex. For s physical this reads 

ReF{s,t) = D{s,t)^ -P.P. r ds'^4^ + - r A{s,t) =ImF{s,t). (2.1) 

JaM^ s' - S TT J4M^ S' -U 

(P.P. denotes Cauchy's principal part of the integral). 

Actually, and because, in some cases, the A[s,t) grow linearly with s, (2.1) is divergent. This 
is repaired by subtractions; that is to say, by writing the Cauchy representation not for F itself, but for 
F{s,t)/{s — si) where si is a convenient subtraction point, usually taken to coincide with a threshold. 
This introduces a function of t in the equations (the value of F{s,t) at s = si); we leave it to the reader 
to rewrite our equations with the appropriate subtraction incorporated 

Let us separate out the high energy contribution, s > Sh (we will fix Sh later) to (2.1). We then 

have 

Dis, t)^- P.P. r ds' + - r ds' + V{s, t; Sh) (2.2a) 

TT i4M2 S' - S TV JiMl S' - U 

and 

V{s,t;sh) = - ds — ^ '-+ ds — ^ '-; (2.2b) 



TT 



D(^\s,t) = 2^- ^{21 + l)PKcos^) . (2.3b) 



I Sh " " ■■' Sh 

we are assuming s < Sh- Both D and the A may be written in terms of the same set of phase shifts by 
expanding tliem,"^ for fixed s channel isospin /, as 

2„l/2 °o 1 

cot^Sl'{s) + l 

One of the factors 2 above occurs because of the identity of pions; we work in the limit of exact isospin 
invariance. 

These equations provide constraints for the phase shifts provided one knows (or has a reliable 
model) for the high energy term, V{s, t; Sh). They enforce analyticity and s ^ u crossing symmetry. 

2.2. The Roy equations 

Eqs. (2.1) to (2,3) look rather cumbersome. Roy'^' remarked that they appear simpler if we project them 
into partial waves, integrating over physical {t < 0) values of the cosine of the scattering angle: one finds 
the Roy equations 

T!(n/ =E/ <^s'Ku,{s,s') + ms;s,). (2.4) 

cot^ 5Y'{s) + l JiMl cot^ 4 ^ (s') + 1 

Here the kernels Ku' are known and the Vi are the projections of V . 

Eq. (2.4) is valid in the simplified case we are considering here, i.e., without subtractions. If we 
had subtractions, the fixed t dispersion relations would acquire an extra term, a function g{€) (the value 
of F{si,t) at the subtraction point). This may be eliminated, using crossing symmetry, in favour of the S 
wave scattering lengths. Eq. (2.4) would be modified accordingly. 



^ We are actually simplifying in that (2.2) should take into account the different isospin structure of s and u 
channels, which the reader may find in e.g. the text of Martin, Morgan and Shaw.''^' 
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Let US rewrite the Roy equations in the form 

C = m,V) (2.5) 

where ^ = {I™/(}^o stands for the set of imaginary parts of the partial waves, for s < Sh, and is the 
functional that follows from (2.4). We can define a mapping, 

e = mV) (2.6) 

and then the solution of the Roy equations is a fixed point of <P. 

The relations (2.5) are highly nonlinear integral and matrix equations. Solutions are known to 
exist in some favorable cases; in fact, Atkinson!*^! proved, even before the advent of Roy's equations, that, 
for any arbitrary V{s, t; Sh) such that it is sufficiently smooth and decreasing at infinity, one can obtain, 
by iterating (2.6), a solution not only of the Roy equations, but of the full Mandelstam representation, 
and compatible with inelastic unitarity for all s as well. Therefore, the solutions to the Roy equations are 
ambiguous in an unknown function, and the matter of what is an acceptable V becomes crucial. This is 
particularly so because fulfillment of the Roy equations does not guarantee full analyticity and crossing; 
and it may happen that a given solution of the Roy equations is incompatible with other sum rules (as is 
the case for the CGL solution). 

3. The scattering amplitude 



At low energy, say s^/'^ < 0.82 GeV, the inelasticity in nn scattering is known experimentally to be 
negligible; it is for these energies that the Roy equations (2.5) are to be solved. To do so we need 
as input the function V or, equivalently, the imaginary part of the scattering amplitude for energies 
gi/2 > Q 82 GeV. In fact, for the Roy equations we need Im F(s,t) for s physical and t physical, t < 0. 
However, for other applications, we will require ImF(s, t) up to the edge of the Martin-Lehmann ellipse,^ 
t < 4M^; our discussion will also cover this case. We now proceed with a discussion of the different waves 
and energy regions. 

3.1. The S and P waves for E below 0.82 GeV 

Because we want to test the solution of CGL for the nn S matrix, we consider now the solution to the 
Roy equations, incorporating chiral perturbation theory to two loops, given there. The low energy SO, S2, 
P waves are written by these authors as 

tan4^)(s) = A;2Vl-4M2/s{4 + Bfk^M^ + Clk^/M^ + Dlk^/Ml) ^MlzJlL^ (3.1a) 

k = a/ s / 4 — , and the values of the parameters, as given by CGL, Eq. (17.2), are 

Ag = 0.220, B§ = 0.268, C§ = -0.0139, L>g = -0.139/10^, sqo = 36.77 

yl^ = -0.444/10, = -0.857/10, = -0.221/10^, = -0.129/10^, so2 = -21.62 (3.1b) 
^1 = 0.379/10, Si = 0.140/10^^, Ci = -0.673/10", = 0.163/10^, si= 30.72 M^. 

These are the values of the phase shifts that we will use up to the energy E = 0.82 GeV. 

To test dispersion relations, either in the form of the Olsson relation or the Proissart-Gribov 

representation, we need also the values of the S, P waves at intermediate energies (0.82 < E < 1.42) and 
the values of the D, F waves below 1.42 GeV, that we take from experiment; higher waves are presumably 
negligible. Moreover, we require ImF(s,t) for s^l'^ > 1.42 GeV. This last we will obtain from Regge 
theory in Subsect. 3.4; we now turn to the intermediate energy regions. 

Before doing so, however, we want to emphasize that, in the present paper, we do not deal with the 
matter of the consistency of the fits for the S, P waves between 0.82 and 1.42 GeV that we will give in next 

For analyticity properties of tttt scattering see for example ref. 7. 
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subsection, or of those for the D, F waves. These fits [Eqs. (3.2) to (3.11) below] are merely a convenient 
way to summarize the experimental data; our results would change very little if we had instead used a 
spline interpolation for the experimental phase shifts. We will discuss this further in Subsect. 4.4.1, where 
we will show that the discrepancy will remain essentially unchanged provided we demand a resemblance 
to the data (allowing for a large uncertainty) in this intermediate region. 

3.2. The S, P waves between 0.82 and 1.42 GeV 



For the SO wave in the region between 0.82 GeV and KK threshold we use the parametrization, obtained 
by fitting experimental data^ (as in ref. 6), 



cot 6^^>{s) = 



Ml 



Ml 



2k .s-iM2 M2 



Bo 



sl/^ = 2Mk; xVd-O-f- = 11-1/(19 - 4). 
M^ = 806 ±21, i3o = 21.91 ± 0.62, Bi = 20.29 ± 1.55, B2 = 22.53 ± 3.48; 

4°' =(0.226 ±0.015) M-\ 



(3.2) 



The solution depends on the value of 6q^^ (Mj^) we impose in the fit. In (3.2) we took that following 

from the more recent measurements of if27r decay, PI dj^^ [M^] = 41.5±3° ; another possibility is to average 

this with the older determination, thus imposing the value = 43.3 ± 2.3° ; this we will discuss 

in Subsect. 5.2. 

Solution (3.2) is, up to 5^2 ^ o.84 GeV, similar to the CGL one, (3.1); see Fig. 3.1. We will use 
the CGL solution up to 0.82 GeV, slightly above their nominal maximum range, s^/2 = 0.8O GeV, and 
(3.2) between 0.82 and 0.96 GoV. 

For the S2 wave between 0.82 GeV and 1.42 GeV we use the phase shift obtained by fitting 
experimental data and including the requirement 0"^^ — 0.044 ± 0.003 (this last follows from the 
analysis of CGL): 



(2),. r 

cot 6\ ' (s) = — ^—^ <Bo + Bi —= , > 

° ^ ^ Ik s-1zl \ ^ Vs + ^fs^^s J 



(3 3) 

s^^ = 1.45 GeV; xVd-o.f. = 16.1/(18 - 2). 

So = - 115 ±4, Si = -106 ±3, ^2 = 139.57 MeV [fixed]. 

This actually corresponds to ag — —0.0457 ± 0.0074. One can allow to vary by 8 MeV, and still be 
within 1 a of the minimum, but we will not do so here. 

Then we have the P wave between 0.82 GeV and 1.0 GeV. Here we fit the pion form factor, 
including e+e" and r decay data. There are now two possibilities: the first one is 

s^^ = 1.05 GeV; xVd.o.f. = 1.3. (3.4) 
Mp = 772.3 ± 0.5 MeV, Bq = 1.060 ± 0.005, Bi = 0.24 ± 0.06. 
[0.82 GeV< <1.0 GeV]. 



® The (slight) differences with some of the parameters in ref. 6 occur because now we are using iW^ = = 
139.57 MeV instead of the average pion mass, 138 MeV, and we are also essentially eliminating from the fit the 
data for energies above 0.96 MeV. The change in the x^/d.o.f. of the 7 = S wave corrects an error there. We 
send to this reference for details on the fitting procedure. 



-4- 



-ON THE PRECISION OF CHIRAL-DISPERSIVE CALCULATIONS OF TTTT SCATTERING- 




FlGURE 3.1. The 1 = 0, S-wave phase shifts corresponding to (3.2) (con- 
tinuous Hne) and Colangelo, Gasser and Leutwyler, ref. 2 (dashed hne). 



In particular, for the low energy parameters, this gives 

fli = (40.6 ± 1.3) X 10"=^ M-^, h = (4.18 ± 0.43) x 10"^ M^^ 

This result is obtained from the fit to the pion form factor, with only statistical experimental errors taken 
into account, performed in ref. 10. If we also take systematic normalization errors into account, (3.4) is 
replaced by 

,1/2 



cot 5i(s) = 1^ (M^ - s) |so + 



and now 



s^^ = 1.05 GeV; xV^-o-f. = 1.1. (3.5) 
Mp = 773.5 ± 0.85 MeV, Bq = 1.071 ± 0.007, Bi = 0.18 ± 0.05. 
[0.82 GeV < < ^ q GeV] 

ai = (38.6 ± 1.2) X 10"^ M'^, bi = (4.47 ± 0.29) x 10"^ M-^ 
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Figure 3.2. Fits to the 7 = 0, 5*- wave phase shift and inelasticity from 
960 to 1350 MeV. Also shown are the data points from solution 1 of Pro- 
topopescu et al.'^^' (black dots) and some data of Grayer et al.'^^' (open 
circles) . 



We will consider both possibilities, but the calculations of dispersive and Froissart-Gribov integrals will 
be made with (3.5), for definiteness. If using (3.4) the differences would be minute.^ 

We next turn to the SO, P waves in the higher energy regions, but still E < 1A2 GeV. For 
the SO wave between KK threshold, 0.992 GeV, and 1.42 GeV, we use a semiempirical formula that fits 
reasonably well the existing data,'^^'^^! from s^/^ > 0.96 GeV to 1.50 GeV: 

,.(0), , {s - M^){M] - s)\k2\ V3-4M^ 
cot6y{s)=co j^j^, ; h = ^^ 

'^^l- ci-r7^+C2— , (3.7) 



sl/2 s J S 

[0.992 < < 1 42 GeV] cq = 1.36 ± 0.05, ci = 6.7 ± 0.15, C2 = -17.6 ± 0.7, 
Mk = 496 MeV, = 0.802 GeV, Mf = 1.32 GeV, M' = 1.5 GeV . 

Note that, for inelastic scattering, we define our parameters so that, in general, 

51 '{s) ^ 

in the elastic region, r}f\s) = 1. The fit to the data following from (3.7) is shown in Fig. 3.2. 



The fact that the errors for ai, bi are smaller when using (3.5) than when using (3.4), which at first sight appears 
counterintuitive, can be understood as follows. The errors in the parameters Bo, Mp in (3.5) are larger than those 
in (3.4) -as suggested by intuition. The error in Bi, however, is smaller, and it is his quantity that influences 
most 61 (the error in oi stays essentially constant). Including systematic errors makes the determinations of the 
pion form factor in the timelike and spacelike regions more compatible one with another, and this allows a more 
precise determination of low energy parameters. 
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Figure 3.3. Fits to the / = 0, D-wave phase shift. Also shown are 
the data points from solution 1 of Protopopescu et al. '^^1 (black dots) and 
some data of Estabrooks and Martin'^^' (open circles). 



Finally, for the P wave between 1 GeV and 1.42 GeV, we use an empirical formula, obtained 
adding a resonance (with mass 1.45 GeV) to a nonresonant background: 

1 Mlr^ 

Im h (.s) = ,,,,,,,, , + BR- ' 



1 + [A + l.lfc2/.sl/2]2 (..,. _ m2 )2 + M^r-^ ' 

[1.0 < < 1.42 GeV] Mp> = 1.45 GeV, T = 0.31 GeV, A = 2.6 ± 0.2; BR = 0.25 ± 0.05. 

(3.8) 

Note that the effect of the p(1450) is very small, as will be clear in our various evaluations below. 
3.3. The D, F waves below 1.42 GeV 

We take these waves as given (from threshold to 1.42 GeV) by the fits of ref. 6, with inelasticity added for 
the DO wave and, for the F wave, including also the tail of the ps resonance. Moreover, we have required 
(for compatibility with the CGL analysis) that the corresponding scattering lengths agree within errors 
with those given in CGL; that is to say, we include the CGL values, weighted with their errors, in the fits 
for D2, F (for the DO wave it is not necessary, as there are enough precise experimental data). For the DO 
wave we thus write 

cot^f (s) = ^ (M/, - s)M^ jso + Si^^^S} ; = 1.430 GeV; 

M/, = 1270 MeV , Bq = 23.7 ± 0.7, Bi = 22.9 ± 2.7. (3.9) 

77 = 1-2x0.15: 



2[fc/A:(M2j]io 



l+[/cA(M2J]20- 
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Figure 3.4. 


Fits to the / 


= 2, 


D-wave 


phase 


shift. Also shown are the 



data points of Losty et al. (open circles) and from solution A of Hoogland 
et al. (black dots), refs. 13. 



The inelasticity on the /2 is taken from the Particle Data Tables. (3.9) corresponds to Og^^ = 
(15 ± 3.5) X 10-^ against CGL's value (17.5 ± 0.3) x 10"^ M'^. For the D2 wave/ 

So^^ = 1.43 GeV; Bq = (2.33 ± 0.17) x 10^ Bi = (-0.39 ± 0.75) x 10^, Z\ = 90 ± 11 MeV . 

(3.10) 

Now = (1.6 ±0.4) X 10-^M^5 [CGL's value: (1.7 ±0.13) x 10"* Af-^]. 
The D phases are depicted in Figs. 3.3, 3.4. 

Finally, for the F wave we write a background plus a Breit-Wigner. The background is obtained 
fitting low energy, the resonance is the ps with its properties taken from the Particle Data Tables: 

Im/3(s)=— — ^ BR- . 



cot ^3(5) 



l + cot^Js \k{Mp,)J (s-MD^ + Min' 



2k^ 

Mp,,_ = 1.69 GeV, T = 0.161 GeV, BR = 0.24; 
Bq = (1.07 ± 0.03) X 10^ Bi = (1.35 ± 0.03) x 10^ 
Here 03 = (7.0 ± 0.8) x 10"^ M'^; the value reported in CGL is (5.6 ± 0.2) x 10"^ M'^. 



(3.10) corrects a mistalce in the corresponding wave in ref. 6. 
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Pb, 



Figure 3.5. Cut Pomeron ladder ex- 
changed between the partons pA and ps 
in hadrons A, B. The emitted gluons 
will materialize into a shower of particles. 
The cross section is proportional to the 
square of the cut ladder. 



3.4. High energy: the Regge picture 



As we will discuss in Sect. 5, the experimental phase shift analyses become unreliable as soon as the 
inelasticity is large; for tttt scattering, this occurs at and above E ~ 1.4 GeV. Fortunately, Regge pole 
theory provides an input for high energy scattering; we will now briefly describe those of its features 
that are of interest to us. Before starting with the details, however, it is perhaps worth while to remark 
that Regge theory is as much part of QCD as, say, chiral perturbation theory; in fact, Regge theory is 
probably of more general validity than QCD. By using Regge formulas we are thus not introducing extra 
assumptions. The only debatable point is when is Regge theory applicable; QCD only specifies s ^ A^, 
s ^ \t\. Fortunately, factorization allows us to relate tttt to ttN and NN cross sections. From this, and the 
fact that Regge formulas and experimental cross sections for tttt scattering agree (within errors) around 
s^^^ — 1.4 GeV, as shown in Figs. 3.6-3.8 below, we will conclude that Regge formulas are applicable at 
and above these energies; specifically, we will use them above E = 1.42 GeV. We now turn to a brief 
discussion of the details. 

Consider the collision of two hadrons, A + B A + B. According to Regge theory, the high 
energy scattering amplitude, at fixed t and large s, is governed by the exchange of complex, composite 
objects (known as Regge poles) related to the resonances that couple to the t channel. Thus, for isospin 1 
in the t channel, high energy scattering is dominated by the exchange of a "Reggeized" p resonance. If no 
quantum number is exchanged, we say that the corresponding Regge pole is the vacuum, or Pomeranchuk 
Regge pole; this name is often shortened to Pomeron. In a QCD picture, the Pomeron (for example) will 
be associated with the exchange of a gluon ladder between two partons in particles A, B (Fig. 3.5). The 
corresponding formalism has been developed by Gribov, Lipatov and other Russian physicists in the 1970s, 
and is related to the so-called Altarelli-Parisi, or DGLAP mechanism in deep inelastic scattering. I^^l 
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Figure 3.6. The average cross 
section ^[2a^o^+ +fT^o^o], which 
is pure 7* = 0, arbitrarily nor- 
malized. Broken line: experimen- 
tal cross section. Note that the 
bump here, as the larger bumps 
in Figs. 3.7, 3.8, is due to the 
coincidence of two resonances, 
/o(1270), /2(1370), mostly elastic, 
around s^/^ ~ 1.3 GeV. Thick 
gray line: Regge formula (3.17a). 
The thickness of the line covers 
the error in the theoretical value 
of the Regge residue. 



One of the useful properties of Regge theory is factorization; I^^l it can be proved from general 
properties of Regge theory.^ Factorization states that, for example, the imaginary part of the scattering 
amplitude FA+B^A+B{s,t) can be written as a product 

lmFA+B^A+B{s,t) fAit)fBms/s)''^^'K (3.12) 

Here s is a constant, usually taken to be 1 GeV^ (we will do so here); the functions /a, /b depend on 
the corresponding particles (if we had external currents, also on their virtuality), but the power (s/s)""^*) 
is universal and depends only on the quantum numbers exchanged in channel t. The exponent afi{t) is 
the Regge trajectory associated to the quantum numbers in channel t and, for small t, may be considered 
linear: 

afl(t) ~ afl(O) + a'flt. (3.13) 

For the p and Pomeron pole, fits to high energy processes give 

ap{0) = 0.52 ± 0.02, a'^ = 1.01 GeV"^ 

ap(0)=l, a'p = 0.11 ±0.03 GeV 

The Regge parameters taken here are essentially those in the global fit la of Rarita et al.t^^l; for ap(0), 
however, we take the value 0.52 ± 0.02 which is more consistent with recent determinations based on deep 
inelastic scattering.^ The results depend very little on this. 



(3.14) 



In potential theory the proof can be made mathematically rigorous; in relativistic theory, it follows from extended 
unitarity or, in QCD, in the DGLAP formalism, as is intuitively obvious from Fig. 3.3. 

Fits to deep inelastic scattering processes, and references to previous literature, may be found in the book by 
FJY in ref. 14 
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Let US consider the imaginary part of the tt-ZV or NN scattering ampUtudes (here by NN we also 
understand NN). We have, 

ImF(f^(.,t) ^ [fl^'\t)]\s/sr-^'\ lraFi%\s,t) c fi'^\t)fl^'\t){s/sr-'^'l (3.15a) 
For If — 1, R — p; for If — 0, R — P (the Pomeron). Therefore, using factorization, we find 

ImFii^\s,t) ~ [/(^*)(f)]' (s/s)"«(*^ (3.15b) 
The functions f\^''\t) depend exponentially on t for small t and may be written, approximately, as^° 

/(^'=°)(t)=a,(P)e^*, /f*=i)(f)=a,(p)i±^ [(1 + 1.48)e''* - 1.48] ; 6 = (2.38 ± 0.20) GeV-^ 

(3.16) 

The exponent b appears to be the same for rho, Pomeron and P' , within errors. '^^1 
Prom (3.15) we can deduce the relations among the cross sections 

CTrAT^all (^NN^all 

and from these relations one can obtain the parameters (Jt^ in (3.16) in terms of the known -kN and NN 
cross sections. Using this, we can write explicit formulas for tttt scattering with exchange of isospin /( = 
in the t channel: 

ImF(^*=°)(s,i) I 1 +0.24W- I a,(P)e''*(s/s)«^(°)+"'^*, (3.17a) 

t fixed I V A' J 

and we have added empirically the subleading contribution, proportional to y/s/s, of the so-called P' pole 
(associated with the /2 resonance) that is necessary at the lowest energy range (see Pig. 3.6). Por It = 1, 



\mFi':Z'l{s,t) ImFW(s,t)+ImF(^'^)(s,i), 

- + ap{t) 
. + ap(0) 



ImFW(s,0 =a^(p)l±^ [(1 + 1.48)e''* - 1.48] (s/s)"''(°)+<*, (3.17b) 

/s \ 1/2 

ImF(^'^)(s,t) =(0.4±0.1) f - j ImF('')(s,t). 

We have added a background (Bk) contribution to the isospin 1 amplitude; this should be considered purely 
empirical and is adjusted so that the asymptotic formula joins smoothly the experimental amplitude at 
low energy, within errors; see Pig. 3.7. 

Prom (3.16) and the known cross sections for nN , NN scattering we have 

cr^(P) = 3.0±0.3; (T^(/9) =0.84 ±0.10 (3.17c) 

where the errors are obtained by considering the dispersion of the values of the parameters in ref. 15, 
and increasing the result by 50%, which should cover amply the uncertainty on the point where one joins 
experimental and asymptotic formulas (that here we have taken to be 1.42 GeV) as well as errors in the 
parameters we have taken fixed. 

It is important to note that the Regge parameters in the fit of Rarita et al. are obtained by global 
fit to ttN , NN and NN data for small momentum transfer and for cm. kinetic energies in the region 
between 1 GeV and 6 GeV (approximately) , which is the region of interest for us here as the contribution 
to the various integrals above this energy is negligible. The results of Rarita et al. are still the best 
available as indeed there are essentially no new data in that energy range. We have, on the other hand. 



^'^Consistency requires a more complicated form for the residue functions fl^^\t)\ see refs. 4, 15. For the small 
values of t in which we are interested, our expressions are sufficiently accurate. 
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Figure 3.7. The cross sec- 
tion a^^*~^\ for isospin 1 in the 
t channel, arbitrarily normalized. 
The dotted line is experiment; 
the short-dashed line the rho ex- 
change Regge theory. The long- 
dashed line is obtained by adding 
to this the p(1450) contribution. 
Finally, the thick gray line in- 
cludes also the background Regge 
piece. The thickness of this line 
is equal to the error due only to 
the rho Rcggc residue (the to- 
tal error for the full theoretical 
formula used in the text, that in- 
cludes errors due to p(1450) and 
the background Regge piece, is 
some 20% larger). 



verified that the cross sections are compatible with the corresponding values as given in the more recent 
editions of the Particle Data Tables. I^^l 

We will treat the errors in the various Regge parameters as uncorrelated. In fact, the leading 
Regge amplitudes (Pomeron and rho) are uncorrelated; there is some correlation with, respectively, the 
It = 2 exchange (see below) and the Bk piece for It = 1 exchange, because they have been fixed by 
fitting the sum to the pion cross sections. Since this only affects subleading pieces this would only have 
a minute influence in the results (in fact, they would slightly decrease the overall error due to the Regge 
contributions), and anyway the variations are substantially smaller than the 50% increase in the errors of 
the Regge residues with which we have made our evaluations. 

For each individual process n^n'^, 7r°7r°, we have to incorporate the amplitude for exchange of 
isospin /( = 2 in the t channel, which would be due to double rho exchange. This cannot be obtained from 
factorization, since nN or NN do not contain such amplitude. We use an empirical formula, 

ImF(^*=2)(s,i) = C2e- 

and we have obtained the constant C2 by fitting the difference between the experimental 7r°7r° and Tr^Tr"*" 
total cross sections at s^/^ = 1.42 GeV, and the Pomeron plus P' values; see Fig. 3.8. 

The dependence of our results on lm.F^^^\ lm.F^^*~-^^ is very slight (for the second, with the 
exception of the 62^^). 

We now add a few words on the matter of when one may apply formulas like (3.17,18). From the 
QCD, DGLAP version of the Pomeron, we expect the following pattern to occur: in the region \t\ <C s, 
s ^ (with A ~ 0.3 GeV the QCD parameter) the ladder exchange mechanism will start to dominate 
the collision A + B. We then will have the onset of the Regge regime with, at the same time, a large 
increase of inelasticity and a smoothing of the total cross section according to the behaviour (3.17). 

For ttN, NN scattering this occurs as soon as one is beyond the region of elastic resonances; in 
fact (as can be seen in the cross section summaries in the Particle Data Tables) as soon as the kinetic 
energy or lab momentum is above 1 to 1.2 GeV. For tttt we thus expect the Regge description to be valid 
for the corresponding energies, that is to say, for s^/^ ^ 1.4 GeV. Indeed, around s^/^ ~ 1.4 GeV it is 
still possible to calculate the nn scattering amplitudes reliably from experimental phase shifts and indeed 




C2 = 0.8 ± 0.2, 



(3.18) 
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3 



2 - 




J I I I I J 

1.1 1.2 1.3 1.4 1.5 (GeV) 



Figure 3.8. The cross sections cr(7r°7r+) (dashed hne), o-(7r"7r°) 
(dotted hne), and the Pomeron plus P' (continuous line). The thick 
gray bands are obtained including the 7 = 2 exchange contributions. 
Their thickness corresponds only to the error of the Pomeron piece. 

they agree, within a 10%, with the Regge expressions in the tt^tt cases; see Figs. 3.6, 3.8. Moreover, the 
experimental inelasticity for tttt around 1.4 GeV, ^ 20%, also agrees with the value of the inelasticity 
measured at the same energies for ttN or NN scattering. 

For the It = 1 amplitudo, and because it is a difference between large amplitudes, the influence 
of resonances may be expected to extend to higher energies. Indeed, we see in Fig. 3.7 that agreement 
between experiment and the Regge expression (within errors) around 1.4 GeV requires adding the resonance 
p(1450), as in Eq. (3.8). We will do so in our calculations. Thus, for all tttt amplitudes we will assume the 
Regge formula (eventually adding the jo(1450) contribution) to be valid for s^^^ > 1.42 GeV. 

As is clear from this minireview, the reliability of the Regge calculation of high energy pion-pion 
scattering cannot go beyond an accuracy of ^ 10%, even for small t. The deviations off simple Regge 
behaviour are expected to be much larger for large \t\, because the counting rules of QCD imply a totally 
different behaviour for fixed t/s. This is one of the problems involved in using e.g. the Roy equations 
that re([uirc integration up to — t ^ .s ~ 1.7 GeV^, where the Regge picture fails completely (we expect 
instead the Brodsky-Farrar behaviour, CTfixed cose ~ s^)- However, for forward dispersion relations or the 
Proissart-Gribov representation we will work only for f = or i = 4M^ for which the largest variation, 
that of e''*, is still small, since b x (t = 4M^) ~ 0.19. So we expect no large error due to departure off 
linearity^^ for the exponent in f-^*\t) or for the Regge trajectories, anit). 

4. Olsson's sum rule and the Froissart-Gribov calculation of ai, bi, aj^^ 63^^ 



"^"'in the case of the rho trajectory, exact linearity would imply (1:^(0) = l/2(Mp — M^) ~ 0.87 GeV not far from 
the value 1.01 GeV~^ that the actual fits give, and which we have used here. 



- 13 - 



-J. R. PELAEZ AND F. J. YNDURAIN- 



4.1. The Olsson sum rule 

The Olsson sum rule is simply a forward dispersion relation for the amplitude F'-^'^^^ with isospin 1 in 
the t channel, evaluated at threshold. Expressing F(^'^^)(4M^, 0) in terms of the scattering lengths, this 

(0) (2) f°° , ImF(^*=i)(s,0) , , 

- 5a^o^ = Do, Do = / ds ); 4.1 

In terms of isospin in the s channel, 

F(^'=i)(s,t) = ^F^'^=°\s,t) + ^F^'^=^\s,t) - |F(^==2)(s,t); (4.2) 
the F^^") are normalized by 

ImF(^=)(s,i) = V(2Z + l)P;(cos9)Im//^^'(s). 

TT/C ^ 

Substituting in the right hand side above the S, P phases of CGL up to 0.82 GeV, the phases (as given 
in the parametrizations of Subsects. 3.2,3) for the same at intermediate energies (0.82 < E < 1A2 GeV), 
the D, F phases from (3.9-11), the tail of the p(1450) resonance between 1.42 and 1.6 GeV, and the p plus 
background Regge parameters of Subsect. 3.4 we find, for 2oq°^ — ba^^^ in units of M^, 



CGL, direct CGL, dispersive 

0.400 ± 0.007 [CGL S, P, s^^ < o.82 GeV] 

0.146 ± 0.004 [Rest, s^/^ < 1.42 GeV (incL, D, F below 0.82 GeV)] 

0.073 ±0.010 [Regge, p s^''^ > 1.42 GeV] 

0.010 ± 0.003 [Regge, Bk; sV2 > 1.42 GeV] (4-3) 

0.005 ± 0.001 p(1450), 1.42 < s^a < i.g GeV 

0.663 ± 0.007 0.635 ± 0.014 [Total, disp.] 



By "direct" we mean the value of the corresponding quantity (in our case, as given 

in CGL. By "dispersive" we understand that we have used the dispersive formula, Dq in (4.1), to calculate 
the same quantity. The "Rest" are the contributions of the D, F waves below 1.42 GeV, plus the S, P 
waves between 0.82 and 1.42 GeV. Of this "Rest", the largest contribution comes from the DO and P 
waves. 

The error in the CGL S, P piece below 0.82 GeV we obtain by varying the A, B, C, D parameters 
in (3.1) according to the formulas given by ACGL (Appendix). It is almost identical to the error given 
for the whole of the "direct" quantity itself. We will discuss in some detail the discrepancy between the 
"direct" and "dispersive" determinations of this quantity 2aQ°' — 5a^^ as the situation for the other as and 
6s to be considered below will be very similar. 

The reason the analysis of the discrepancy is not straightforward is that both determinations are 
strongly correlated, as they both depend on the same parameters. The "direct" determination is obtained 
from the parameters A, B, C, D in CGL (as given in Eq. (3.1) here), which describe in particular the S 
waves. So we should really write 

fOfl^o) - ^„ (2)1 "direct" 

The integrals in the dispersive determination contain the contributions of the S, P waves up to 0.82 GeV, 
which are given by the same A, B, C, D, so one also has 

(0) . (2)/'dispersive, S, P CGL" 
[zuq ouq \a,b,c,d 

Now, it is clear that if we try to change the parameters A, B, C, D in the "direct" determination to, for 

example, decrease the value of a^^ to bring this closer to the dispersive value, the same change in the 
A, B, C, D will affect the integral over the SO wave in [2ao°^ — ba'^^]^^^^^^^'^'^'^^ , making this smaller 
and therefore nullifying to a large extent the improvement. 
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What one has to do to solve this problem is to consider the difference A = ^af^^ — 5oq'^^ — Dq and 
vary here the parameters A, B, C, D. Then we find the value 

A = 0.027 ±0.011 : 

that is to say, a 2.5 a discrepancy. 

This procedure will also be followed for the Froissart-Gribov sum rules, where the correlation in 
the CGL analysis is transmitted in part by the common chiral perturbation theory parameters li. (We 
will discuss more about errors in Subsects. 4.4 and 5.1). 



4.2. The Froissart-Gribov representation: ai, bi 



By projecting the dispersion relation (2.1) (or a derivative with respect to t of it) over the Ith partial wave 
in the t channel, at i = 4M^, one finds the Froissart-Gribov representation 

V^rjl + l) lmF{s, 4M2) 



4M^r(? + 3/2) " s'+i ' 

V^m + l) p /4ImF^„^,(s,4M2) {I + l)lraF{sAM^)' 



' 2M^r{l + S/2) J^M2 \ (s-4M2)s'+i s'+2 

ImF^^^g = {d/dcos9s) ImF. For amplitudes with fixed isospin in the t channel, an extra factor 2 (due to 
identity of particles) has to be added to the left hand side; so we have, for example, 

o„(/=i)_ V^r{l + 1) p ^ JmF(^^='HsAM^) 



AM^r{i + 3/2) y4M2 

With the same type of calculation as for the Olsson sum rule, and with the same definitions, we 
now find, in units of 10~^ x M~^, 



ai, CGL, direct CGL, Froissart - Gribov TY (St.) TY (St.+Sys.) 

18.5 ± 0.2 [CGL S, P, < o.82 GeV] 
9.1 ± 0.3 [Rest, §1/2 < QeV] 
8.1 ±1.1 [Regge, p] 

1.0 ±0.3 [Regge, Bk] (4-5) 
0.3 ±0.1 p(1450) 

37.9 ±0.5 37.1 ±1.3 [Total, Froissart-Gribov.] 40.6 ± 1.4 38.6 ±1.2. 

Here, and for hi , we profit from the existence of an independent determination of the P wave parameters, 
using the pion form factor data both in the timelike and in the spacelike regions, ^^^^ denoted by TY. From 
this we have chosen two values: from the fit taking into account only the statistical errors in the various 
data sets (St.), as in Eq. (3.4) here; or taking also into account the systematic normalization errors (St. 
+ Sys.), as in Eq. (3.5). 

The distance between the direct evaluation and the one with the Froissart-Gribov calculation is 
now 0.6(7, and there is also acceptable overlap with the TY (St.+Sys.) figure. 
For the quantity bi we have, in units of 10~^ x M~^, 
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61, CGL, direct CGL, Proissart - Gribov TY (St.) TY (St. + Sys.) 

-0.92 ± 0.05 [CGL S, P, s^/^ < 0.82 GeV] 
1.02 ± 0.04 [Rest, s^/^ < GeV] 
5.33 ± 0.86 [Regge, p] 

0.55 ±0.16 [Regge, Bk] (4-6) 
0.01 ±0.0 p(1450) 

5.67 ±0.13 5.99 ±0.88 [Total, Proissart-Gribov.] 4.18 ±0.43 4.47 ±0.29. 

Here the Regge contribution is particularly important because the lower energy pieces cancel almost 

completely. The numbers labeled TY, as before, refer to what one obtains from the fit to the pion form 
factor. We remark that this last is a very robust determination in that it is obtained by fitting some 210 
points from several independent experiments, is independent of high energy assumptions and it covers 
spacelike as well as timelike momenta: thus, the values of the threshold parameters are obtained by 
interpolation, notoriously more stable than extrapolations. 

There is no inconsistency between the "direct" and Proissart-Gribov numbers for the CGL calcu- 
lation, but they are both too large by almost 4(7 compared to even the more favorable value, TY(St.+Sys.), 
following from the pion form factor. 

4.3. The Froissart— Gribov representation: a2^\ 62^'; 7 = 0, 2 

We first calculate the two combinations of scattering lengths ao+ = §[02°^ "^2'^^] ^^'^ <^oo — §[0^2'^^ 
They correspond to the s— channel amplitudes 

The only important difference with the cases in the previous subsection is that the dominant high energy 
part is given now by the Pomeranchuk trajectory (instead of the rho) and its importance is small because 
the integrals converge faster. We find, in units of 10~^ x M~^, 

ao+5 CGL, direct CGL, Proissart — Gribov 

8.43 ± 0.09 [CGL S, P, < 0.82 GeV] 
1.84 ± 0.05 [Rest, s^/^ < QeV] 
0.68 ± 0.07 [Regge, It = 0] 

-0.06 ± 0.02 [Regge, It = 2] (4-7) 
0.04 ±0.01 [yo(1450) 

10.53 ±0.10 10.94 ±0.13 [Total, Proissart-Gribov.] 

In finding the error of the "direct" value, (10.53±0.10) x 10~^ M~^, it is important to take into account the 

strong correlations of the errors of the 02'^'', a'^P . To do this, we use Eq. (14.4) in ACGL to calculate directly 
the quantity ao+. The difference between the "direct" and Proissart-Gribov values, with correlations taken 
into account, as we did in the case of the Olsson sum rule, is now 



0.38 ± 0.09, 



so that the discrepancy reaches the 4 a level. 

In the same units, 10~^ x M~^, we have 
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aoo, CGL, direct CGL, Proissart — Gribov 

11.73 ± 0.32 [CGL S, P, s^/^ < 0.82 GeV] 
1.91 ± 0.04 [Rest, s^/^ < 1.42 GeV] 
0.68 ± 0.07 [Regge, It = 0] 
0.12 ±0.04 [Regge, It = 2] 

13.94 ±0.32 14.44 ±0.33 [Total, Proissart-Gribov] ; 



we have also taken into account the correlations a la ACGL to evaluate the error of the "direct" number. 
The difference between "direct" and P.-G. values for CGL are, with correlations taken into account, of 

0.49 ± 0.09, 

i.e., a 5a discrepancy. 

Pinally, we present the results for 5o+ — |[&2^^ — ^2^''] and 6oo — §[^2*^^ + 2&2'^^], both in units of 
10-^ X M-^: 



6o+5 CGL, direct CGL, Proissart — Gribov 

-0.331 ± 0.015 [CGL S, P, s^/^ < 0.82 GeV] 
0.04 ± 0.00 [Rest, .s^/^ < 1.42 GeV] 
0.12 ±0.02 [Regge, It =0] 
-0.05 ± 0.02 [Regge, It = 2] 

-0.189 ±0.016 -0.233 ±0.036 [Total, Proissart-Gribov.] 



The contribution of the resonance /?(1450) is now negligible. Por the difference between the direct and 
Proissart-Gribov result we have 

0.044 ± 0.026, 
that is to say, almost a 2(T discrepancy. Por boo, 



boo, CGL, direct CGL, Proissart — Gribov 

-6.90 ± 0.22 [CGL S, s^/^ < 0.82 GeV] 
0.07 ± 0.01 [Rest, s^/^ < 1.42 GeV] 
0.12 ±0.02 [Regge, It = 0] 
0.10 ±0.05 [Regge, /t = 2] 

-6.72 ± 0.22 -6.62 ± 0.23 [Total, Proissart-Gribov.] 



Por 6oo the direct result and the one following from the Proissart-Gribov representation differ by 2 a: 

0.10 ±0.05. 

However, one cannot take this or the discrepancy for bo+ as seriously as in the previous cases. This is so 
because of the large (relative) size of the contribution of the /t = 2 exchange piece, proportional to the 
derivative with respect to t of an expression we have obtained purely empirically by fitting at t = 0. 
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4.4. How significant cire the discrepancies? 



In the present subsection we investigate whether the inconsistencies we have found can be eliminated (or 
to what extent they can be made less severe) by altering the non-CGL part of the dispersive, or Froissart- 
Gribov calculations. We will do so in two steps. First, we will consider what happens if we alter the pieces 
labeled "Rest" in (4.3) to (4.10); then we will address the question of what can be done at high energy 
(sV2 > 1.42 GeV). 

4.4.1. The region between 0.82 and 1.42 GeV 

We start with the first question, that we discuss in detail for the Olsson sum rule since the results for the 
Froissart-Gribov calculations are very similar. We then consider the following set of drastic modifications 
of our calculations: For the SO wave, and 0.82 < E < 0.992 GeV we may replace (3.2) by the CGL 
parametrization, (3.1). For the SO wave and 0.992 < E < 1.42 GeV, where it is poorly known, we allow 

5q'^ to vary between the two extreme values tt and 37r/2. For the S2 wave, we multiply by 3 the errors 
given in (3.3). For the P wave, and 1 < E < 1.42 GeV, we change the elasticity of the p(1450) resonance 
by 50% (up and down). For the DO wave, that supplies the more important contribution to "Rest", we 
consider the effect of taking the /2(1270) resonance to be purely elastic, or 30% inelastic. The remaining 
contribiitions to "Rest" are so small that we need not worry about them. 

The alterations just discussed are rather extreme; nevertheless, their effects are of no relevance. 
They produce the following extra errors (we give the central value of each term as well): 

SO, 0.82 < < 0.992 GeV : 0.026i^:E]o6 
SO, 0.992 < < 1 42 QgV : 0.018l[j;[]?^ 
S2, 0.82 < < 1.42 GeV : -0.022 ± 0.004 

P, 1.0 < 5^/2 < 1 42 GeV : 0.024 ± 0.005 
DO, < 1 42 GeV : 0.055 ± 0.001. 

Including these increased errors we get that, for the Olsson sum rule, the result for the "Rest" changes 
according to 

"Rest" : 0.145 ± 0.004 ^ 0.145lg;g?^, 
and, for the whole dispersive result, we now get 

Total : 0.631 ± 0.013 ^ 0.63lig gl^ , 

i.e., practically no change at all in the upper error bar. 

4.4.2. The high energy region, s^/'^ > 1.42 GeV 

Once we have verified that the inconsistencies between the CGL direct and dispersive calculations of low 
energy parameters cannot be due to the contributions of the intermediate energy region, we turn to the 
high energy (s^/^ > 1.42) piece. Then, we relax the condition of factorization for the p and Pomeron Regge 
residues (but we do not change the others) . We treat them now as free parameters, describing an effective 
scattering amplitude, to see under which conditions one can reconcile the direct and Froissart-Gribov (or 
dispersive) evaluations for the scattering lengths and effective range, in the CGL-like analysis. Starting 
with the isospin 1 case, we thus write 

lmF^^\s,t) AcT.(p)l±^ [(1 + 1.48)e''* - 1.48] {s/sr^('^+<\ 

f fixed + Oip\^) 

that is to say, we modulate the p amplitude in (3.17b) by the constant A. We then fix a^{p) — 0.85, and 
treat A as a free parameter. We then find that overlap between the direct and dispersive determinations 
for the quantity 2ao°^ — ba'g^ involved in the Olsson sum rule would require A = 1.4, which is well outside 
expectations and, moreover, this spoils the overlap for oi, bi, which become inconsistent at the 2 to 2.5 cr 
level. 
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For the the situation is even more transparent. Consider for example the quantity ooo, 
Eq. (4.8). Integrating only to 0.82, with the CGL phases, we find 11.73, which is the bulk of the result. 
Even if the errors of what we call "Rest" were underestimated by a factor 3, and this "Rest" would be 
1.79 (instead of 1.91), adding it one would get at least 13.52 ± 0.33 for the contribution below 1.42 GeV. 
The direct result, with the CGL values of the af , is 13.94. To get agreement, one would require the high 
energy, E > 1.42 (Regge) estimate to be wrong by a factor 2, very difficult to believe. And it would be no 
good: the same Pomeron that contributes to aoo contributes to ao+ and to the 6o+) ^oo- The disagreement 
would be shifted to the ^^oo, which would then be wrong by about 4(7, and ao+ would still be wrong 
by almost 2a. As for the proverbial square peg in the round hole, trying to fit a corner only makes others 
sick out more sharply. 

5 Discussion of the ACGL and CGL analyses 

5.1. Possible cause of the distortion of the CGL solution 

In this section we try to ascertain the reasons for the troubles that seem to afflict the CGL analysis. This 
is particularly important because, althoiigli ACGL or CGL did not verify the Froissart-Gribov relations, 
they did check relations similar to the Olsson sum rule. It follows that the reasons for the discrepancies 
must be due to the high energy input. Here you have two regions: between 0.82 and 1.42 GeV (more or 
less) the inelasticity is low, and, as we have shown, one can trust the experimental phase shifts. Even if 
they have systematic errors, these will likely not be large and they will just produce a slight fluctuation 
of the solution of the Roy equations, as we have shown explicitly in Subsect. 4.4.1 that it occurs for our 
evaluations. 

The difhcult region, however, is for s^^"^ above 1.42 GeV. Between 1.42 and 2 GeV, CGL pre- 
sumably use the phase shifts of ref. 11 and, above 2 GeV, a Regge-type formula. We start the discussion 
with the region 1.42 < s^^'^ < 2 GeV. Here inelasticity is very high, and the phase shifts and inelasticity 
parameters cannot be determined reliably, at the level of accuracy required. Of course, you can always 
give numbers that fit the experimentally observed moments in peripheral two-pion production; but so will 
other, in some cases very different values of 5s. and r/s. In the energy region 1.4 < 5^/2 < 2 GeV, the phase 
shifts and inelasticities all stem from a single set of experiments and are likely to disagree with reality 
by much more than their nominal errors. In fact, this can be seen to occur for the S wave even at lower 
energy: as soon as the KK channel opens, the Cern-Munich phase shiftsl^^l disagree violently with the 
Berkeleyl^'^] ones. This emphasizes the dangers of relying on a single experiment for the phase shifts, as 
one has to do already for s^/^ > 1.2 GeV. 

It is not difficult to see how different phases may give similar results, for the elastic cross section. 
For example, consider the elastic wn cross section, in the P wave: in both cases (Cern-Munich and Particle 
Data Tables results) it is small. In the Cern-Munich one, because sin^ Si is small; in the other because r] 
is small. Unfortunately, the imaginary parts of the inelastic amplitudes are very different; contrary to the 
Cern-Munich results, in the PDG case it would be large, at least around the resonances, because of the 
contribution of the inelastic channels. The converse (i.e., overestimate of the total cross section) may, of 
course, also happen. In fact, the cases mentioned are just examples of an ambiguity (over and above that 
due to experimental errors) proved to exist quite generally in ref. 16, and which is likely to be large as 
soon as you have important inelastic channels open. 

Now, CGL (following Penningtonl"*^^! ) take the Cern-Munich phase shifts, that probably contain 
large and unknown systematic errors, and impose sum rules [e.g., the sum rules (B.6,7), (C.2) in ACGL], 
following from low energy crossing symmetry, to fix the Regge parameters at energies E > 2 GeV. Not 

""^^This unreliability is reflected, for example, in the Particle Data Tables (e.g., the edition of ref. 15), where no 
number is given for the branching ratios of tttt resonances with masses at or above 1.2 GeV (with the exception 
of the p3(1690)) and even the SO phase around the /o(980) has a dubious status: this last due to the ambiguity 
caused by the inelastic KK channel. In fact, the resonances that appear in tttt production are not the same that 
one finds in e~''e~, r decay or J/ip decay, and the inelasticities in both cEises are also quite different. 
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surprisingly, CGL (and Pennington^^) get irrealistic Regge parameters (as realized by CGL themselves); 
for example, ACGL and CGL get a Pomeron with a width of the diffraction peak which is s-independent, 
and twice the standard value (at low s), and a residue much smaller than what factorization implies. In 
fact, we will show in the Appendix explicit calculations of two sum rules (in particular of the sum rule 
(B.7), one of the crossing sum rules that Pennington and ACGL use) which are perfectly satisfied by a 
standard Regge amplitude, with factorization for the rho and Pomeron trajectories, provided one uses 
Regge asymptotics from s^^"^ = L42 GeV. 

According to CGL this deviation from conventional Reggeistics is not important because the 
influence of the high energy region (s^^^ > 1.42 GeV) into their low energy (s^/'^ < 0.82 GeV) phase shifts 
is very slight. However, and as we have shown in the present paper, inconsistencies show up as soon as one 
considers sum rules -like the Proissart-Gribov sum rules- that are sensitive to the high energy behaviour 
of the amplitudes. 

From the previous analysis it thus follows that CGL start, in the Roy equations, from a V with 
incorrect Regge behaviour and dubious phase shifts above 1.42 GeV. Let us call this F(Wrong R). CGL 
run this through the Roy equations (2.6) and find a solution, ^(WrongR). Now, this solution is not 
horrendous because experimental low energy data, chiral perturbation theory and crossing sum rules force 
you to have the errors in Regge parameters and cross sections compensating, to a certain extent, in what 
regards their low energy effects. Indeed, the independence on the low energy partial waves on the high 
energy amplitudes used is approximately true, for the ACGL results, where the mismatch that occurs if 
using the correct Regge asymptotics stays below the 2 a level. However, for the CGL results, use of chiral 
perturbation theory (with neglect of higher order corrections) has the dual effect of highly correlating the 
various low energy parameters and excessively decreasing the errors. Thus, for example, the value for the 
quantity a+o that follows from the Proissart-Gribov representation, 10.94 ± 0.13 (in units of 10~^M~^) 
is displaced 4a from the value following directly from the parameters of CGL, 10.53 ± 0.10. Now, a+o 
is directly related to the chiral constant I2, 0+0 — [h — 27/20]/7207r"^/^M^. Hence a variation of a+o 
implies a corresponding variation of I2 , or of higher chiral corrections, that destabilizes all the quantities 
that depend on it in a chiral perturbative analysis; in particular, the low energy S and P waves. As we 
have shown in the present paper, inconsistencies show up in the CGL scattering amplitude (with standard 
Regge parameters) as soon as one considers sum rules that, like the Proissart-Gribov or Olsson ones, are 
sensitive to the high energy behaviour. What the inconsistencies found in the previous section show is 
that the distortion is several times larger than the nominal CGL error bars. 



5.2. A tentative alternate solution 



In support of the idea that the effects discussed in the previous subsection are indeed the cause of the 
mismatches in the CGL 5* matrix, we have calculated the Olsson sum rule and the quantities ai, &i using 
now, for s^^^ < 0.82 GeV, the results of the fit, wave by wave, reported in ref. 6, Sect. 7.6. Por the wave SO 
we take now the fit obtained imposing the value ^q*^^ (M^) — 43.3±2.3° , and with only three parameters; 
we then have 

Bo = 21.04, Bi = 6.62, = 782 ± 24 MeV; (5.4a) 

/-F=T^; 4°^ = (0.230 ±0.010); ^^(^if) = 41-0° ± 2.1° . 
Q.o.i . iy — o 



^^This should not be taken as a criticism of the work of Pennington; at that time the data, very poor, did indeed 
suggest possible deviations of Regge behaviour for tttt. On the other hand, the fact that QCD also implies 

factorization was of course unknown. 

^''This fit is actually a refinement of that of Eq. (7.6.2) in ref 6; more details about this will be presented in a 
separate publication. 
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FlGURE 5.1. The / = 0, S'-wave phase shifts corresponding to (5.4) (continuous line) 
and Colangelo, Gasser and Leutwyler, ref. 2 (dashed hne). Some experimental points are also 
shown. 



The errors of the Bi are strongly correlated; uncorrelated errors are obtained if replacing the Bi by the 
parameters x, y with 

Bo=y-x; Bi = 6.62 - 2.59a;. (5.4b) 

Then, 

t/ = 20.04 ± 0.75, x = 0±2.4. (5.4c) 

The solution is shown, compared to the CGL phase, in Fig. 5.1. We then integrate with (5.4) up to 
E = 0.82 GeV and with (3.2) from 0.82s to KK threshold. For S2, P we take the same fits as before, 
specifically, eqs. (3.3), (3.5). 

We find the following results, in units of M^: 

Olsson direct dispersive . 

0.691 ± 0.042 0.659 ± 0.020 ^ 

[here "direct" means that we take the values following from the fits in Eqs. (3.3), (5.4)]. Moreover, and 
also in units of M„, 



ai direct, (TY, St.) direct, (TY, St. + Sys.) Froissart - Gribov 
(40.6 ± 1.4) X 10-3 (38.6 ± 1.2) x lO'^ (37.9 ± 1.4) x lO'^ 

and 
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direct, (TY, St.) direct, (TY, St. + Sys.) Proissart - Gribov 
(4.18 ± 0.43) X 10-3 (4.47 ± 0.29) x lO^^ (5.69 ± 0.96) x 10"' 



The tag "direct" now refers to the values of ref. 10, with only statistical errors (St.) or including also 
systematic errors (St. + Sys.). Thus, we find agreement at the la level in all three cases; for ai, bi, with 
the TY (St.+Sys.) solution. With the same parameters we find, for the D waves, and with the help of the 
Proissart-Gribov representation, the values 

ao+ = (10.60 ± 0.17) X lO"'' aoo = (14.99 ± 0.68) x lO"'' M^^ (5.8a) 

and 

bo+ = (-0.170 ± 0.083) X lO"'' M"^, 6oo = (-6.91 ± 0.47) x 10"^ M"^. (5.8b) 

This is compatible with what we found for the Oj^' by direct fit to the experimental data in Sect. 3.3 
within the rather large errors of these last values. 

The large error, and the separation in the central values in the Olsson sum rule, Eq. (5.5), is due 

to the fact that the data do not fix with sufficient accuracy the Og scattering length, which provides most 
of the error in the "direct" number. In fact, as is known, one can use the Olsson sum rule to refine the 
parameters of the S2 wave; if we do so, fixing all other parameters to their central values (within errors) 
and include the Olsson sum rule in the fit to the S2 wave we find 

cot5^^)(.)=£^^ |j?o + i?,^-^- 



1/2 o (5-9) 

4' = 1.45 GeV: x /d.o.f. = 17.2/(19 - 2). 

B„ = - 118 ± 2.5, Di = -105 ± 2.5, Z2 = 139.57 MeV [fixed]. 
Then one has a^^^ = -0.0428 ± 0.0022 and (5.5) becomes 

Olsson direct dispersive 

0.671 ±0.023 0.663 ±0.018. ^ ' 

The rest of the relations (5.6-8) improve slightly, and the D wave scattering lengths also change a little: 

ai = 38.0 ± 1.2 X 10-3 M'^, hi = 5.64 ± 0.96 x lO'^ M'^; 

ao+ = (10.51 ± 0.15) X 10-^ M'^, aoo = (14.89 ± 0.65) x 10"^ M'^. 

It should be noted that the error here for ao+ is at the edge of the region of credibility, as indeed it is of 
the order of magnitude of electromagnetic corrections which the analysis does not take into account. This 
value of ao+ implies, at one loop level, a very precise value for the chiral perturbation theory parametert^^l 
I2 of 

I2 = 5.97 ± 0.07 

Of course the agreement in (5.5,6,7,10) is not enough to guarantee that the new solution is consis- 
tent; to prove that, one would have to check the whole set of dispersion relations and crossing constraints, 
something that will be the subject of a separate paper. But it clearly suggests that the CGL solution fails 
to pass the tests because it is distorted. This can also be inferred by comparing the CGL solution for 
the S2 wave with (3.3) as in Fig. (5.2), where we show the CGL and (3.3) together. While both fit the 
data below 0.82 GeV [expression (3.3) gives actually a slightly better fit even there], the distortion of the 
CGL solution above that energy is suggestive. A similar pattern is found in Figs. 3.1, 5.1. This very much 
suggests that the CGL fit is a forced fit, biased by a reflection of a faulty high energy scattering amplitude. 
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FlGURE 5.2. The / = 2, 


iS-wave phase shifts corresponding to (3.3) 



(continuous hne) and Colangelo, Gasser and Leutwyler, ref. 2 (dashed 
hne) . Also shown are the data points of losty et al. (open circles) and 
from solution A of Hoogland et al. (black dots), refs. 13. 

6. Summary and conclusions 

We have checked a number of tests of the low energy (s^/^ < 0.82 GeV) SO, S2 and P wave phase shifts 
given in ref. 2 by Colangelo, Gasser and Leutwyler, based on two loop chiral perturbation theory plus 
the Roy equations with a certain high energy (s^/^ > 1.42 GeV) input. We have shown that, if we used 
the values for this high energy piece that follow from Regge theory, then the Olsson sum rule and the 
combinations of scattering lengths Oo+ = |[ct2°^ — 012'^^], ^00 = |[a2°'' + 'ia^2^] show mismatch by as much 
as 4 ~ 5 (T. We have discussed in detail why we think that the discrepancy is inherent to the low energy 
(s^/2 < 0.82 GeV) CGL phases. Thus, in Subsect. 4.4.1 we have shown that even rather drastic alterations 
of the middle energy region, 0.82 < s^/"^ < 1.42 do not alter the inconsistencies. 

With respect to the higher energy region (s^/'^ > 1.42 GeV), the situation is such that, if one tries 
to modify the Regge piece to fit the Olsson sum rule (say) then not only the alteration (40 to 100%) is 
much more than what one can reasonably expect, but the lack of consistency is shifted to ai, 61. A similar 
phenomenon -in fact, even more pronounced- occurs with 00+ and aoo- This we discussed in detail in 
Subsect. 4.4.2, where it is clear that the mismatch is due to the low energy CGL input. Moreover, the 
value of the quantity bi remains displaced by 4 a from what one gets from a fit to the pion form factor. 

It should be borne in mind that we are talking here about disagreements at the level of a few 
percent; so, if one is prepared to shift the central values of CGL by iij) to 2 a, and double their errors, the 
inconsistencies disappear. This is what happens, for example in the analysis of ACGL where the errors 
are from 3 to 10 times larger than those in CGL. Nevertheless, at the level of precision claimed by CGL, 
the errors are real. We have argued that they are probably due to an irrealistic high energy (s^/^ > 1.42 
GeV) input, which distorts the low energy phase shifts. In support of this we have shown that a direct fit 
to data, including fully analyticity constraints, for the P, SO, S2 waves (in the case of the last, requiring 
also consistency of the Olsson sum rule to decrease its errors) plus a high energy input given by orthodox 
Regge theory, produces a different set of compatible low energy phase shifts and high energy scattering 
amplitude. This set is formed by the phase shifts given in Eqs. (3.5), (5.4) and (5.9). This set is in fact 
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similar to that of CGL, but is shghtly displaced and its errors are slightly larger; so, for example, the 
quantity 5q°^(M|-) — S^^Mj^), important for kaon decays, changes according to 

Si°\M'^)-5^o\M]^)= 47.7 ±1.5° [CGL] ^ 48.4±2.1° [Our solution]. 

A fact that may be mentioned here is that Descotes et a.'^' have, in a recent article, found a 

solution whose central values differ from that of CGL by almost 2(CGL) standard deviations and in fact 
point in the direction of our tentative alternate solution here. Thus, they have, in units of M^r, 

= 0.228 ± 0.012, = -0.0382 ± 0.0038 [Descotes et al.] 

Their errors are also more like what we have in our alternate solution. Note, however, that whether or 
not the alternate solutions turn out to be consistent has nothing to do with the consistency of the CGL 
solution: this last fails independently of the failure or success of the novel one(s). 

Analyticity determines the real part of the scattering amplitude in terms of its imaginary part. 
However, to get the real part you need to know the imaginary part up to infinity. Now, if the imaginary 
part is wrong at high energy and yet the dispersion relation (or Roy equations) are satisfied, it necessarily 
follows that one must have made a compensating error in the low energy imaginary part. In other words: 
you have fallen into a spurious solution. The fact that the solution is spurious should be manifest as soon 
as one devises a test that gives a difFerent weight to high and low energy pieces: this is exactly what we 
do in our paper, for the CGL solution, with the help of the Proissart-Gribov representations. 

Appendix 



In this Appendix we briefly discuss (and prove the failure of) the reason for the unorthodox Reggeistics 

chosen in ACGL, CGL, in as much as it has a bearing on our subject matter here. These authors, following 
Pennington, [^^] set up crossing sum rules [Eqs. (B.7), (C.2) in ACGL], which relate high and low energy, 
and conclude that they are satisfied only if, in particular, the Pomeron residue is about 1/3 of the value 
implied by factorization and the rho residue is about a 40% larger. 

Contrarily to the conclusion of ACGL, however, we will show by explicit calculation of two typical 
sum rules that, if one assumes orthodox Regge behaviour from s^^^ > 1.42 GeV, the low energy phase 
shifts are perfectly compatible with the value of the Rcggc residues implied by factorization. This will cinch 
the proof that, as discussed in Subsect. 5.1, the Reggeistics of ACGL are very likely due to compensation 
of the unrealistic phase shifts used between 1.42 < s^/^ < 2 GeV. 

Specifically, we will consider a sum rule dominated at high energy by the Pomeron, viz., the sum 
rule (B.7) in ACGL; since it is independent of the S and P waves, it constitutes a new, stringent test of 
the Regge structure. We will also consider another sum rule, which is dominated by the rho trajectory. 

The first sum rule may be written as 



f ^4 



Here F'(^) (s, t) = dF^^^ (s, t) /d cos 6, and the index / refers to isospin in the s channel. 
We will separate J into a low energy and a high energy piece: 

j-Sh /-oo 

J = J\.e. + Jh.e.\ J\.e.= j ds . . . , Jh.e. = ds . . . . 

■JAMl Jsh 

The low energy piece, J\.e., only contains contributions of waves D and higher. Since these waves are only 
known with (relatively) large errors, it is (generally speaking) very dangerous to draw conclusions about 



^^These errors are particularly large, and uncertain, above 1.3 GeV, where inelasticity begins to be important. For 
example, the error on the DO wave contribution to Ji.e. due to a 50% change in the inelasticity of the /2(1270) 
resonance is as large as the nominal error due to only the errors in the parameters in (3.9). 
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the high energy integral, Jh.e. , from the experimental value of the low energy piece, Ji.e.- Nevertheless, 
we will show that, if we choose Sh = 1.42^ GeV^, then we find perfect consistency, within errors. In 
this calculation we will first neglect the contributions of exchange of / = 2 and of the backgroimd to rho 
exchange, both of dubious status and substantially smaller than the Pomeron and rho exchange pieces, 
but we keep the P'. Using the parametrizations of Sect. 3.3 for the D, F waves we find, in units of M~^, 

J\.e.{B waves) = 1.222 x 10~^, Ji.e.(F wave) = -0.076 x 10"'' 

so that, including the errors, 

Ji.o. = (1.15 ± 0.05) X IQ-^ (A. 2) 

For the high energy piece, expanding in amplitudes with definite isospin in the t channel, and with the 
numbers in Sect. 3.4 for the Pomeron and rho contributions, we get 

Jh.e. (Pomeron) = -1.093 x 10"'*, A.^Xp) = 0.034 x 10"^, 

i.e., including errors, 

^h.e. = (-1-06 ± 0.17) X IQ-^. (A.3) 

Thus, we have cancellation between (A. 2) and (A.3), within errors: there is no reason to justify departure 
off the expected Regge behaviour. 

We next comment a little on the P' and on the inclusion of the It =2 contribution. Because the 
high energy part of the sum rule (A.l) is mostly given by the t derivative of the even isospin amplitudes, 
a more precise evaluation than the one carried here would require that we replace the P' contribution 
of (3.17a) by a more accurate formula. Unfortunately, the characteristics of this Regge pole are poorly 
known; see ref. 14. If we take for the the P' trajectory a formula like that of the p, then (A.3) is replaced 
by 

Jh.e. (With corrected P') = (-1.2 ± 0.2) x 10"^. 

Including also the Jt = 2 contribution, as given in (3.18), we would find 

Jh.e. (With corrected P', and including = 2) = (-0.5 ± 0.3) x 10"^. (A.4) 

This still cancels the low energy piece, (A. 2) at the 2 a level. This discrepancy cannot be taken seriously, 
because of the uncertainties in the P' trajectory and because the t slope in formula (3.18) is little more 
than guesswork. 

The second sum rule is obtained by profiting from the threshold behaviour to write an unsubtracted 
forward dispersion relation for the quantity F^^^^^\s, 0)/(s — 4M.^) This gives the relation 

^^^v^ifd^^-f"^^^^^^^ ^d.^^^i^^, (A.5) 



Ml 



which is known at times as the (second) Olsson sum rule. The index / refers to isospin in the s channel 
and C^j""* are the s — u crossing matrix elements. Canceling a\ with the Proissart-Gribov expression for 
this quantity and substituting the C^j""* we find the result 

J Ml J Ml s^{s-4M^y 

(A.6) 

The contributions of the S waves cancel in (A.6), so only the P, D and F waves contribute (as 
usual, we neglect waves G and higher). At high energy, I2 contributes little since the corresponding 
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integral converges rapidly: most of the high energy contribution comes from the first term, dominated by 
rho exchange. We will use units so that = 1 and obtain the following results: 

/(low energy, P wave) = (-2.80 ± 0.3f) x 10"^ 

/(low energy, DO + D2 waves) = (0.56 ± 0.03) x 10"^ 

/(low energy, F wave) = (0.01 ± 0.00) x 10"^ 

/(high energy, p) = (2.41 ± 0.37)) x 10"^ ''^''^'^ 

/(high energy, / = 0) = -(0.17 ± 0.02) x IQ-^ 

/(high energy, / = 2) = -(0.02 ± 0.01) x IQ-^ . 

By "low energy" we understand the contributions from energies below 1.42 GeV, where we use phase shifts 
and inelasticities to calculate the scattering amplitudes, and "high energy" is above 1.42 GeV, where a 
Regge description is employed. The final result for the sum rule is 

/= (0.016 ±0.37) X IQ-^ 

i.e., complete cancellation of low and high energy contributions. 

The remarkable fulfillment of these sum rules show the incorrectness of the assertions found in 
ACGL, CGL: both for Pomeron and rho, standard Regge behaviour for wn scattering is perfectly consistent 
with crossing symmetry provided one imposes it systematically for energies above 1.42 GeV. 

Note added in proof 



After this article was sent to the publisher, a preprint has appeared [I. Caprini, G. Colangelo, J. Gasser 
and H. Leutwyler, hep-ph/0306122j in which some of the conclusions (but not the calculations) of our 
work are contested. We do not think it necessary to alter our paper on account of the work of Caprini et 
al.; we plan to present a discussion of it in a separate article. 
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